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In this paper new mult i- variable deformed Fokker-Planck (FP) equations are presented. These 
deformed FP equations are associated with the Ruijsenaars-Schneider-van Diejen (RSvD) type sys- 
tems in the same way that the usual one variable FP equation is associated with the one particle 
Schrodinger equation. As the RSvD systems are the "discrete" counterparts of the celebrated exactly 
solvable many-body Calogero-Sutherland-Moser systems, the deformed FP equations presented here 
can be considered as "discrete" deformations of the ordinary multi-variable FP equations. 

PACS numbers: 05.40.-a, 02.30.1k, 03.65.-w , 02.30.Gp 



The Fokker-Planck (FP) equation is one of the most important tools to deal with fluctuation phenomena in various 
kinds of systems [l[ . Recent discovery of anomalous diffusion in fractal and disordered media have prompted new 
developments in the FP theory, the major one being fractional FP equations, where the ordinary spatial derivatives 
are replaced by fractional derivatives 0. Nonlinear FP equations Q and diffusion equation based on g-derivatives 
have also been considered Q. 

In based on the well-known relation between FP equations and Schrodinger equations of quantum mechanics, 
we have proposed new types of deformed FP equations which are associated with the Schrodinger equations of the 
"discrete" quantum mechanics considered in 0, 0, Q ■ The latter is a natural discretization of quantum mechanics 
and its Schrodinger equations are difference instead of differential equations. The cigenfunctions of some exactly 
solvable "discrete" quantum mechanics include various deformations of the classical orthogonal polynomials (the 
Hcrmite, Laguerre and Jacobi polynomials), namely, those belonging to the family of Askey-scheme of hypergeometric 
orthogonal polynomials, and the Askey- Wilson polynomial in g-analysis [9|, [l0|, ■ In fact, these deformed orthogonal 
polynomials arise in the problem Q of describing the equilibrium positions of Ruijsenaars-Schneider-van Diejen 
(RSvD) type systems [HI, which are the "discrete" counterparts of the celebrated exactly solvable multi-particle 
Calogero-Sutherland-Moser (CSM) systems (HI, just as the Hermite, Laguerre and Jacobi polynomials describe the 
equilibrium positions of the CSM systems !14]. In this sense, the proposed deformed FP equations Q can be considered 
as the "discrete" deformation of the usual FP equation. 

In this paper, we would like to generalize the results in [5( for the one-particle cases to the many-body cases, and 
obtain the many-variable deformed FP equations associated with the RSvD systems. We first briefly review in Sect. II 
the connection between the ordinary multi-variable FP equation and the many-body Schrodinger equation. Sect. Ill 
summarizes the basics of the CSM systems associated with rational/trigonometric classical root systems. In Sect. IV 
we review the RSvD systems and derive the corresponding FP equations for the two rational types. The trigonometric 
type RSvD systems and the corresponding FP equations are discussed in Sect. V. Sect. VI concludes the paper. 



II. CONNECTION BETWEEN FOKKER-PLANCK AND SCHRODINGER EQUATIONS 

The general form of FP equation of the probability density P(x, t) with n variables x = t {x\, x%, . . . ,x n ) is [1] 



The functions Dj(x) and Djk(x) in the FP operator Lpp are, respectively, the drift vector and the diffusion matrix 
(we consider only time-independent case). The drift vector represents the external forces acting on the particles, while 
the diffusion matrix accounts for the effects of fluctuation. 

It is well known that in one-dimension, the FP equation is closely related to the Schrodinger equation. The 
correspondence between the two equations can be shown by transforming the FP equation into the corresponding 



I. INTRODUCTION 
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Schrodinger equation or vice versa [f|, based on a similarity transformation. The similarity transformation is 
related to the ground state of the Schrodinger equation, or equivalently, to the stationary state of the FP equation. 
One thus expects that similar connection is possible between a FP equation with multiple variables and a one- 
dimensional many-body Schrodinger equation. It turns out that such connection is possible only for special forms of 
Dj and Djk |l| . Below we shall briefly review such connection. 

Instead of transforming an assumed form of FP equation to a Schrodinger equation, as is usually done in the 
literature (see e.g. [H), we think it is more direct to see what form of FP equation will emerge by directly transforming 
a known many-body Schrodinger equation. So let us consider a quantum mechanical Hamiltonian of n particles with 
equal masses (for simplicity, we adopt the unit system in which K, the mass m of the particles, and the 'fictitious' 
speed of light c in the RSvD systems are such that h — 2m = mc = 1) 

H = -E^+nx)- (2) 

i—l 1 

Here V(x) is the potential term which include external fields and inter-particle forces. By adjusting the additive 
constant in V(x), the energy of the ground state 4>o can always be made zero, i.e. HtpQ = 0. Under these assumptions, 
the Hamiltonian H is completely determined by its ground state wave function <f> (x). By the well-known theorem of 
quantum mechanics, 4>q{x) has no node and can be chosen real. Hence one can parametrize 4>q by a real prepotential 
W(x): 

<j> (x)=e- W(x) - (3) 
Then the condition H(j)o = enables one to express V{x) in terms of W(x) as well: 

3=1 



Since only derivatives of W(x) appear in V(x), W(x) is defined only up to an additive constant Wq. We choose 
the constant Wo in such a way as to normalize <fio(x) properly, J (j>o(x) 2 d n x — 1. For simplicity of presentation, 
we consider the cases in which the ground state wave functions are square integrable, that is the corresponding FP 
operators have the normalizable stationary distributions. 

The Hamiltonian H can be expressed as a sum of factorized forms, H = A^Aj, where 



dxj dxj 
d dW(x) 



dxj dxj 



Since each term A^Aj in H is positive semi- definite, the condition Aj(f>o — must hold in order to have Htfio = 0. 
This latter condition is a sufficient and necessary condition for Htfio = 0. 

Now we define an operator from the H and 4>q by the similarity transformation 

L FP = -foHto 1 , (6) 

which guarantees the non-positivity of the eigenvalues of Lfp- Using cf>Q in (J3]), one obtains 



jfp 



v g ( d(2W(x)) \ v 
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(7) 



From (jTJ) , it is seen that Lpp is just the corresponding FP operator with drift vector Dj = —2dW/dxj and a diffusion 
matrix proportional to the unit matrix: Djk = 5jk (the constant factor here equals unity which is a result of the choice 
of unit in H). This establishes the required forms of Dj and Djk so that the multi- variable FP equation is related 
to a one-dimensional many-body Schrodinger equation. We see that Dj is in fact defined, as in the one variable FP 
equation, by a drift potential $(x) = 2W{x): Dj = —d$/dxj. 

Both H and Lpp are determined by 0o- The relationship between their eigenfunctions is the same as that for 
the single particle case, except that special care is needed to ensure the reality of the eigenfunctions belonging to 
degenerate eigenvalues. The eigenfunction P m {x) of Lpp corresponding to the eigenvalue — A m is related to the (real 
and normalized) eigenfunction <j) m of H corresponding to A m by P m (x) = 4>o(x)(f> m (x). Here m is, in general, a 
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certain multi-index. The stationary distribution is Po — <Po = ex P( — 2W), which is obviously non-negative, and is 
the zero mode of Lpp: LppP — 0. Any positive definite initial probability density P(x,0) can be expanded as 
Cm^m(^)i with constant coefficients c m 

/OO 
4> m (x)(<f>o 1 (x)P(x,0))<Px. (8) 
-oo 

Then at any later time t, the solution of the FP equation is P(x, t) = cj>o(x) J2 m c m'/ , m(i) exp(— X m t). 



III. CALOGERO-SUTHERLAND-MOSER SYSTEMS 



One class of multi- variable FP equations having the special form discussed in the previous section which has been 
widely discussed in the literature is that related to the celebrated exactly solvable many-body quantum system, 
namely, the CSM models. Such FP equation has found applications in, e.g. random matrix theory [lH ] and others 

M- ' 

Here we consider four types of CSM systems associated with the classical root systems: two rational and two 
trigonometric types. As mentioned earlier, we consider only the cases with explicitly known ground state wave 
functions 4>a(x), which are square integrable. This excludes the other well-known cases of the hyperbolic and elliptic 
potentials as well as the rational cases without the harmonic confining potential. The potential Vcs(x) of the CSM 
systems can be written in terms of the prepotential W(x) as in The explicit forms of the four types of potentials 
Vcs(x) and the corresponding prepotentials W(x) are as follows: 



(i) rational A n _ 



n n , _ -I s 

K*(*)=5>?+ E (I _ x 2 -n(l+g(n-l)), Oa) 

j^k 

1 " 

VK ( X )=2^^ 2 " 5 ^ log\xj-x k \+W , (9b) 

3=1 l<j<k<n 



(ii) rational BC n : 



7=1 V 3 ' j,k = 

j^k 



-2n{g s + \+9u{n-l)), (10a) 
W ( x ) = ?E, X ] - 9M \ lo s\ x 3 ~ x k\ +log|zj +x k \) -gs^2\og\xj\ + Wo, (10b) 



2 

3=1 l<j<k<n 3=1 



(hi) trigonometric A r , 



ir , . v-^ 1 n q n(n 2 - 1) , 

W{x) = -9 log|sinf(^-x fe )|+Wo, (Hb) 

l<j<k<n 
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(iv) trigonometric BC n : 



, r , s 7T 2 f (9s + gL)(gs + 9l - 1) . g L {9L-l 

Vcs {x) = » § I — ^ — + 

Qm(9m - 1) 3m (9m - 1) 



— V , 

i2 fli, V sin 2 £ - x k ) sin 2 £ + x k ) 



2 



\ 2 i 9mi2 



yf[(9s + 2g L +g M (n 1))' + -f% 2 - 1)), (12a) 



W(x) = - g<ivfflog|sm f(a;j -a;fe)| +log|sinf (iCj + £fe)|) 

l<j<fc<n 
n 

- ^ (g 8 log|sin fx, | + g L log|sin f 2^ |) + VF . (12b) 

Here Wo is the constant term necessary for the normalization of the ground state wave function (f>o(x). The 
constant terms in Vcs(%) are the consequences of the expression Q in terms of the prepotential. In the above 
formulas <7,<?s,<7m and g^ are non-negative coupling constants, and L is the circumference of the circle in which x 3 - 
live for the trigonometric cases. In the rational potential cases, (i) and (ii), the coefficient of the harmonic confining 
potential a; 2 is usually written as ^mui 2 . Here we choose 2m = muj = 1 for simplicity. The rational D n model is 
obtained from the rational BC n model by putting gs = in pOa[) and (|10b[) . The rational B rll C n and BC n models 
are all equivalent. The trigonometric D n model is obtained from the trigonometric BC n model by putting gg = g^ = 
in (I12aj) and (|12bjl . whereas the trigonometric B n and C n models are obtained by gL = and gs = 0, respectively. 

If we apply the transformation © to the Hamiltonian of the rational A n _i CSM system, the resulted FP equation 



" d 2 » Q 

.7 = 1 .7=1 J 



9Y,-—\- (13) 



V 



k=i x i Xk 



Multi-variable FP equation of this type was first employed by Dyson as a Brownian motion model of certain ran- 
dom matrix ensemble [151 ] . FP equations corresponding to the other three types of CSM systems can be obtained 
accordingly. 



IV. RUIJSENAARS-SCHNEIDER-VAN DIEJEN SYSTEMS 



Now we would like to derive new types of multi-variable FP equations corresponding to certain discrete defor- 
mations of the above CSM systems; namely, the Ruijsenaars-Schneider-van Diejen (RSvD) systems with the ratio- 
nal/trigonometric potentials. 

The RSvD system is an integrable deformation of the CSM system. The Hamiltonian has the general form 

H = E (y/viw e ^'\A^) + y/W^yfiM-Wx) - to) • (w) 

j=l V / 

Here dj = d/dxj, and Vj(x) is given by 

" f 1 for A 

Vj{x)=W{Xj) R V{X >- Xk)X \v( Xj+ x k ) tor BCJ (15) 

We use the conventional notation that V*(x) is the complex conjugate function of Vj(x). (For an arbitrary function 
f( x ) = En a nX n , a n € C we define f*(x) = J2 n a n x " '• H ere c * is the complex conjugation of a number c S C.) Since 
the operators e ±ldj cause finite shifts of the wave function in the imaginary direction {e ± ' ldj f{x) = f(x\, ■ ■ ■ ,Xj ± 
i, ■ ■ ■ , x n )), these systems can thus be called "discrete" dynamical systems. 
The basic potential functions v[x) and w(x) are as follows: 
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(i) rational A % 



(ii) rational BC n : 



(iii) trigonometric A n _i 



v{x) 
w(x) 

v(x) 



v(x) = 1 — i—, 
x 

w(x) = (ai + ix^j (a,2 + ix^j , 

1 — % — , 

x 

(ai + ix){d2 + «a;)(a3 + ix)(ct4 + ix) 
2ix(2ix + l) 



sinf (x - ig ) 



w(x) = 1, 



(iv) trigonometric BC n 



v(x) 
w(x) 



sin f (x - ig ) 



sinf(x — ig{) sin f (a; 



(16a) 
(16b) 

(17a) 
(17b) 

(18a) 
(18b) 

(19a) 



*ff2) 



sin fx sm | (£_ |) 

cosf (x - ig[) cos£(x- § - i^) 



(19b) 



cos fx cos J (x — §) 

Here the coupling constants a\ , 0,2, 0,3, 04, g, ffoi 5i , 52, <?i and g 2 are assumed to be non-negative and one of a%, . . . , CI4 
must be greater than 1/2. Apparently, the deformed theories have usually more coupling constants than the original 
ones. 

Note that the Hamiltonian (fTi)) can be expressed as a sum of factorized forms, H = J^j AjAj, with 



A] = JV j (x)e-^-Jvnx)e 



-id* 



(20) 
(21) 



It should be remarked that in the small momentum (the large c) limit, in which the shift operators e ±ldj can be 
approximated by differential operators e ±ldj «li idj — |<3|, the RSvD systems reduce to the CSM systems @, H|- 

We want now to derive new multi-variable FP equations associated with these RSvD systems, following the proce- 
dure described before. Let us note here that while the trigonometric cases can be treated exactly in the manner to 
be discussed below, it is more suitable to be treated as multiplicative shift systems instead of additive ones like the 
other two cases. For this reason, we will discuss the two rational cases first. 

The ground state eigenfunctions 4>o(x) of (TT3|) for the rational cases are 



(i) : </> (x) cx 

(ii) : <j> (x) cx 



r(ai + ixj)T(a 2 + ixj) ■ JJ 



n 

n 



l<j<k<n 



T(g + i(xj - x fc )) 
T(i{x 3 -x k )) 



Ila=i r ( a " + ix j) . "Q tt T{g + i{x 3 +ex k )) 



r(2ix i ) 



l<j<k<n e=± 



1 T(i(xj + ex k )) 



(22a) 
(22b) 



where |/(^)| = y f(z)*f(z) for any complex function /(z). These ground states have zero energy: -ff^o — 0. Since 
_ff consists of a sum of positive semi-definite terms AjAj, the condition of zero ground state energy implies Aj<pQ = 
(j = 1, . . . ,n), or explicitly, 



V 3 I ^i.' 



2'"' 



V) xi, x, + -,•■•, x„ W>o xi, X,- + 



^0 [ Xx , • • • , Xj 2 ' ' ' ' ' 



(23) 



6 



This can be verified explicitly using the forms of V(x) and 4>q for the two cases. 

Now, we form the associated FP operator from (jT5J) and </>o in (|22ap and (|22b[) according to the similarity transfor- 
mation (6]). With the help of (23]), we find that 

n 

lfp = -Y, ( e%9j y j (*) + e ~ t9j v j (*) - v i (*) ~ y; (*)) • (24) 

3=1 

This is the general form of FP operator corresponding to the discrete Hamiltonian H in (33]). One has Lpp<fry = as 
a consequence of H(j) — 0. Thus 4>q is the stationary solution of the respective FP equation. 

In the small momentum limit, or equivalently, the large c limit, all x and other parameters in V{x) are considered 
small quantities 0, Eq. ([24]) then reduces to the corresponding FP equation associated with the CSM system. 
Explicitly, keeping in the operators exp(±i9j) only up to the second order terms in the momentum, we get the FP 
operator 

n / q2 g \ 

L *r = E d^ Re v i( x ) + Q^ 2Im V ^ x ) ■ ( 25 ) 
j=i \ ■> 3 J 

Here in Vj only the lowest order terms in x and other parameters shall be retained. 
As an illustration, let us take the case of rational A n _\ type. The potential is 

Vj(x) = (ai + ixj) (a 2 + ixj) JJ (l - i — j *\ . (26) 

In the above-mentioned limit, we have 

n 1 

i?e Vj -(x) = dicta, /m = (ai + 02) a?j — 0x025 , (27) 



fc=i 



and hence the limiting FP operator is 




L FP = a ia2 E £ + 2 E 7^ I ( ^ + £ H" ~ ff E r-^— I I • (28) 



fc=i 



This is just the corresponding FP operator for the A n -i type CSM system in IT3]) . The presence of the scale factors 
arises from the way the Hamiltonian and the potential of the RSvD systems are parameterized. If we take 01 = a 2 — 2, 
and rescale H — > H/4, then the two formulas are the same. 

V. MULTI- VARIABLE g-SHIFT FOKKER-PLANCK EQUATION 

As mentioned before, the trigonometric cases can be treated more elegantly in terms of multiplicative shift type 
of Hamiltonian 0, @]. Let us define zj = exp(2isixj / L), q = exp(— 2tt/L) and clq — q 9 " ■ For any function f(z) 
(z = t {zi, ■ ■ ■ ,z n )) with real coefficients, we have f(z)* = /(z _1 ). Furthermore, we define Dj = Zj-£p- Then q Dj 

is a multiplicative q-shift operator, i.e. q Dj f(z) — f{z\, ■ ■ ■ ,qzj, ■ ■ ■ ,z n ). With these, the Hamiltonian (fl4|) of the 
trigonometric case with (fT5|) , (|19ap and (|19b|) can be recast into (up to an overall constant factor) 

n 

H = YA ^VjWqV'JVjiz-i) + JV 3 {z-^)q-^JV 3 {z) - V 3 (z) - Vj(z^ 1 ) ) , (29) 



j 



where the potential Vj (z) have the form 

1 - a Zj/z k 



M : w*)=n ~ 7 ' (30a) 

kjtj Zj/Zk 
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and the {a a }'s are related to the parameters in (|19a[) and (|19b[) (see jy]): 

(a x , a 2 , a 3 , a 4 ) = (e" 2 ™ ' L , e -2*U»+i/a)/£ _ e -2^A ^-2^+1/2)/^ 

The form of the Hamiltonian (f29|) is so chosen that it looks similar to that of the rational case (fT4|) by the formal 
replacement e~ ldj — > , e"^ — > g~ Dj . However, it should be emphasized that the potential Vj(z) in corresponds 
to V*(x) in |H]), which is due to the fact that < q < 1 (q = e _27r/i ). 
The ground state <j>o satisfying H(/)q = is 



(hi) : <f>o(z;go,q) « 



n 

l<j<k<n 



{a a zj/z k ;q) c 



(iv) : (f>o(z;{a a },g ,q) cx 



n 



,=1 r£=l(Oa*j;9)oo !<,U<neSi (°0*i**;«)= 



(31a) 
(31b) 



Here we use the standard notation (a; g)^ = H^o(^ — a 1 n )- 

Just like the rational cases, the Hamiltonian is a sum of factorized forms, i.e. H = J^j AjAj, with 



Ah 



q^^Viiz-^-q-^^z), 



The ground state </>o is annihilated by Af. Aj(f>o(z) = (j = 1, . . . ,n). Explicitly this equation reads 



(32) 
(33) 



-i -l 
q 2 Zj , 



,z n 1 )4> Q (zi 1 --- ,q*zj,--- ,z n ) 



= ■■ ■ ,Q *Zjr-- ,z n )(j) (z u -- ■ ,q 'Zj,--,Zn). 

Using H, (j>o, and ([34]) , the similarity transformation ([6]) produces the FP operator 

n 

L FP = -'£ l (q- D >V i (z) + q^Vjiz- 1 ) V^z) - V^)) . 

3=1 



(34) 



(35) 



This is the general discrete g-deformed FP operator corresponding to the Hamiltonian (|29|) . Again, Lpp annihilates 

A2 



VI. SUMMARY 



In this paper we have proposed new types of multi- variable FP equations associated with the RSvD systems which 
are difference instead of differential equations. As the RSvD systems are the integrable "discrete" counterparts of the 
celebrated exactly solvable many-body CMS systems, the deformed FP equations presented here can be considered as 
"discrete" deformations of the ordinary multi-variable FP equations. They are also integrable. In the small momentum 
(large c) limit, our FP equations reduce to the usual FP equations associated with the respective CMS systems. For 
future study, it would be of interest to investigate solutions of these generalized FP equations with potentials V(x) 
and V(z) other than those exactly solvable ones discussed here. 
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